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00 . Abstract 

o ■ 

t— I ■ This paper investigates the late-time behaviour of certain cosmo- 

logical models where oscillations play an essential role. Rigorous re- 
suits are proved on the asymptotics of homogeneous and isotropic 
""^ ■ spacetimes with a linear massive scalar field as source. Various gen- 

q-i eralizations are obtained for nonlinear massive scalar fields, /c-essence 

', models and f(R) gravity. The effect of adding ordinary matter is dis- 

cussed as is the case of nonlinear scalar fields whose potential has a 
■ degenerate zero. 

X' 

1 Introduction 



Scalar fields are important in cosmology as a mechanism for producing models 
with accelerated expansion, both in the very early universe (inflation) and in 
more recent epochs. The mathematical properties of spatially homogeneous 
solutions of the Einstein equations coupled to different types of scalar field 
and ordinary matter have been studied in a number of papers pi], [5], [9], [TO] . 
[IT] . [3] in the case that there is accelerated expansion in the whole future of 
some late time. A corresponding analysis of situations where the accelerated 
expansion is only temporary, being succeeded by decelerated expansion, has 
not yet been carried out. One reason for this is that in contrast to what is 
found for models with continuing accelerated expansion the behaviour of the 
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simplest models, which are homogeneous, isotropic and spatially flat, does 
not carry over to more general homogeneous models. In this paper the very 
simplest case of this type, the massive linear scalar field, is analysed. Then 
a number of other cases showing similar late-time behaviour are discussed. 

Consider solutions of the Einstein equations coupled to a massive scalar 
field which are homogeneous, isotropic and spatially flat. The standard pic- 
ture of the dynamics of these solutions is as follows PQ. (See also the earlier 
papers [H] and |15j.) The solutions exhibit accelerated expansion at in- 
termediate times and eventually enter a phase where they are on average 
decelerated and resemble dust solutions. The main aim of this paper is to 
prove rigorous results about the late-time behaviour. Formulae for this can 
be found in the literature. (See for instance [5].) Here a proof will be pre- 
sented that these formulae do provide asymptotic expansions for all solutions 
in a precise sense. A curious fact, which does not seem to be mentioned in 
the literature, is that it is a consequence of the asymptotic expressions that 
each solution has infinitely many phases of accelerated expansion. The early 
and intermediate time behaviour will not be discussed further here but note 
that some rigorous results on this were obtained in [8]. 

A key reason why the massive scalar field is more difficult to handle than 
types of scalar field which have previously been treated mathematically is the 
following. When deriving asymptotic expansions in this case it is repeatedly 
necessary to estimate integrals with upper limit infinity which are convergent 
but not absolutely convergent. The convergence of the integral, which is 
required for the proof, results from the cancellation of positive and negative 
contributions. In other words, there are oscillations in the solution which 
continue to have a significant effect for all times. 

The paper is organized as follows. In Section 2 asymptotic expansions are 
proved for the basic example, the massive scalar field. This is generalized to 
a large class of nonlinear massive scalar fields in Section 3. It is shown that 
a number of features of the asymptotics are unchanged and the expansion is 
derived up to the point where the first modification occurs. It is discussed 
briefly how these results can be applied to f(R) theories of gravity. Analo- 
gous results for certain fc-essence models are also derived. In all of this the 
only matter present is the scalar field. The question of the incorporation of 
ordinary matter such as a perfect fluid is discussed in Section 4. In Section 
5 it is investigated what happens for a scalar field whose potential has a 
degenerate minimum where it vanishes. Conclusions and possible extensions 
of the results obtained are the subject of the last section. 
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2 The massive scalar field 



This section is concerned with a massive linear scalar field minimally cou- 
pled to the Einstein equations under the assumptions of a spatially flat Fried- 
mann model and the absence of any other matter fields. To simplify the 
computations it is assumed that the scalar field has unit mass. The scale 
factor is denoted by a and the Hubble parameter is given by H = a/ a. The 
equation of motion of the scalar field is + 3if0 + = and the Hamiltonian 
constraint is H 2 = 4^(0 2 + 2 ). The evolution equation for H is H = — 47T0 2 . 
If H is zero at any time t\ then the initial data for the scalar field at time t\ 
vanishes. Thus the scalar field vanishes at all times as does H . This is just 
flat space and so from now on it will be assumed that H never vanishes. It 
must have a constant sign and this sign is chosen to be positive corresponding 
to an expanding model. The equation of motion of can be rewritten as 

+ v / l^(0 2 + 2 ) 1/2 + = O. (1) 

The function H cannot increase and so it must tend to a limit H as 
t — > oo. It will now be shown that H vanishes. Assume that H > 
with the aim of deriving a contradiction. The functions and are both 
bounded as is H. From the equation of motion it follows that the same is 
true of 0. Differentiating the equation repeatedly with repect to time shows 
that higher derivatives of are also bounded. Let {t n } be a sequence tending 
to infinity and let 4> n (t) = <j>(t + t n ) and H n (t) = H(t + t n ). Using the bounds 
already listed, the Arzela-Ascoli theorem [13] can be applied. This implies 
that, after passing to a subsequence, n converges uniformly on compact sets 
to a limit 0oo. Moreover the first and second derivatives of n converge to 
the corresponding derivatives of 0^. The sequence H n converges uniformly 
to the constant value H . With this information it is possible to pass to the 
limit in the equation (TjQ) to obtain 

0oo + 3#o0oo + 0oo = (2) 

From the Hamiltonian constraint it follows that Hq = ^-(0^ + 2 >o )- Differ- 
entiating this relation with respect to t and substituting (|2J) into the result 
shows that 0^ = 0. Then using d2J) again shows that 0^ = and this leads 
to a contradiction to the assumption that H ^ 0. Thus in fact H(t) — * as 
t — > oo. 

Next more detailed asymptotics will be obtained. The pair 

(0/(0 2 + 2 ) 1 / 2 ,0/(0 2 + 2 ) 1 / 2 ) 
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defines a function of t with values in the unit circle. Thus it is possible 
to define a real-valued function 6{t) whose projection to the circle under 
identification modulo 2tc is the given function with values in the circle. The 
function is unique up to shifts of its argument by integer multiples of 2n. 
Let r = (</> 2 + 2 ) 1 / 2 . Then the following equations can be derived: 



r 




-\/l27r 2 sin 2 

- 1 — Vl27rr sin cos 



(3) 
(4) 



Given e > there exists a time t\ such that yY2mr < e for t>t\. It follows 
that on that interval \9 + 1| < e. This means that 

0! - (1 + e)(t - h) < 9(t) <9 l -(l-e)(t- h) (5) 

where Q\ = 9(t\). The evolution equation for r can be rewritten in the form 



r-^t) = r - 1 (t 1 ) + V3n 



(t-t 1 )+ / cos 20' 't l (0') dff 



e(t) 



(6) 



Suppose that e < 1. Let k be the largest integer such that 9\ — kit > 9{t). 
Then 

" cos26'0- 1 (6')de' 



9(t) 



< Trfl 



(7) 



On the other hand 



01-fcTT 



cos26'6-\6')d6' 



—kir 



cos20'[0- 1 (0') + l]d0 / 



< 
< 



&7re(i -e)" 1 



(8) 



Thus the integral on the right hand side of is bounded in modulus by 
'1 — e) _1 (7r + e(l + e)(t — ti)). Putting this information into © shows that 



r-\t) > r" 1 ^) + >/35F{-(l - e)-\ + [1 - e(l + e)(l - e)" 1 ]^ - t x )} (9) 

As a consequence r _1 grows linearly with t and r(£) = 0(t _1 ). The evolution 
equation for then implies that + 1 = 0(t _1 ) and that 



0(t) = - t + 0{\ogt) 



(10) 
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Putting the improved estimate for 9 back into (jH]) shows that 

r(t) = ^- + 0(r 2 \ogt). (11) 



Translating this to the original variables it can be seen that the leading order 
contribution to H is 2/3t and that the scale factor behaves asymptotically 
like t 2 / 3 . Thus the late-time dynamics is similar to that of a model with dust. 
To go further it is necessary to obtain an improved estimate for 9. Note that 

— (r cos 29) = r cos 29 - 2r sin 299 = 2r sin 29 + 0{r 2 ) (12) 

Integrating the equation for 9 and using the above relation gives 

9 = - t + C + 0{r 1 ). (13) 

for a constant C. This constant can be eliminated by choosing an appropriate 
origin for t. Then 9 = — t + 0(t _1 ). 

Consider now the asymptotic behaviour of the scalar field. 

cos(0(t)) = cos t cos(0(t) +t) + sin t sin(0(t) + t) 

= cost + Oir 1 ) (14) 



Hence 



m = ^-+o(t- 2 \o g t) (is) 



Similarly 

<P(t) = ^l + 0(t- 2 \ogt) (16) 
\/3nt 



Putting the expressions for sin 9 and cos 9 into the equation for 9 gives: 
9 = -1 -r 1 sm2t + 0(r 2 \ogt) (17) 
This can now be substituted into ([12]) to get 
d 2 2 

— (rcos26>) = 2rsin26» =t" 2 sin 2 1 cos 2t + —=r 2 sin 2 2t + 0(r 3 \ogt) 

dt V37T V37T 



(18) 

The following lemma is useful for estimating the integral of oscillatory quan- 
tities. It is proved by integrating by parts twice. 
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Lemma 1 Let F be a smooth periodic function with mean zero, / = F' and 
k a positive real number. Then 

f(t)r k dt = F(T)T~ k + C + 0{T- k - 1 ) (19) 

for a constant C . 

Integrating f|T8|) . applying Lemma 1 and noting that due to the informa- 
tion already obtained the integration constant must vanish, it follows that 

3 + 2 cos 2t 



9{t) 



-t 



At 



+ 0(t~ 2 \ogt) 



(20) 



Putting the information now available into the equation for the time deriva- 
tive of r _1 gives 



d 



and after integration 



3tt 1 — cos 2t + sin It 



3 + 2 cos 2t 
2t 



o(r 2 logt) (2i) 



Sir 



sin 2t 



for a constant t\. It follows that 

h 2 



3(t - ti 



1 + 



2 

sin 2t 



O(rMogt) 



2f 



+ 0(r 2 iogt) 



(22) 



(23) 



Note that the constant ti in the leading term in the expression for H could 
be got rid of by a translation in t at the expense of introducing a similar 
constant elsewhere. The expressions for the scalar field can now be improved 
to give 



and 



i 



^7r(t - ti) 



3 + 2 cos 2t 



COS t - MI! I ( ——- ) I ; ^— + 0(r 3 logt) 

(24) 



sin2t 



2f 



3vr(t - 



3 + 2 cos 2t 
it 



sin 2t 



2f 



si i i / - (■■( »s -/- ( - — — ) i + + o (r* log t) 

(25) 

The first of these may be compared with equation (5.45) on p. 240 of [6]. 

Since the scale factor grows slower than linearly the expansion is on av- 
erage decelerated. On the other hand d = ^-((p 2 — 2<fi 2 ) and so there are 
infinitely many intervals on which d > 0. 
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3 Nonlinear scalar fields 



The aim of this section is to investigate to what extent the results for the 
massive linear scalar field can be extended to the case of a potential of the 
form V(4>) = 4> 2 /2 + W(<j>) where W is smooth and W(<j>) = 0((f) 3 ). The 
equation of motion for the scalar field is 

<j> + v / 12^(0 2 + <j) 2 + 2W((f))) 1/2 <j) + + W'(<f>) = (26) 

The evolution equation for H is unchanged. It can be concluded that H 
tends to a limit Hq as t — > oo. To constrain the value of Hq a compactness 
argument can be used as before, leading to the limiting equation 

0oo + 3# o 0oo + 0oo + W'iM = (27) 

with Hi = ^(0^ + + 2W(0oo))- It follows that ^ = 0. Assume now 
that V has no critical points other than the origin in an interval [— <fii,<fii] 
and that the initial value of H is smaller than \[^f<t>i so that <f> remains in 
that interval. Then it follows that (p^ = and this gives a contradiction. 
Thus in fact H = 0. As in the linear case polar coordinates (r, 9) can be 
introduced in the (0, </>)-plane. In these variables the equations read 

r = -v / i^r 2 sin 2 fl(l + 2r- 2 iy(rcosfl)) 1/2 - sin0W(rcos0) (28) 
= -1 - v 7 !^ sin cos 0(1 + 2r- 2 W(rcos9)) 1/2 

-r~ l W'{r cos 9) cos 9 (29) 

Note that 9 + 1 = 0(r) so that it can be concluded as in the case W = 
that 9 grows linearly with t. The equation for r can be rearranged to give 

j^r' 1 ) = Vl2^ sin 2 9(1 + 2r~ 2 W(r cos9)) 1/2 + sm9r- 2 W'(r cos9) (30) 

The last term in this equation looks worrying but note that up to a remainder 
of order r it is equal to hW"'(0) cos 2 9 sin 9. The last expression has mean 
zero and so its integral up to time up t is a bounded function of t, as can be 
seen by changing the variable of integration from t to 9. It can be concluded 
as in the linear case that r _1 grows linearly with t and that r(t) = 0(t _1 ). 
Furthermore 9 + 1 = 0(t _1 ) and (flQj) and (iTTj) hold. Thus H(t) has the same 
leading order behaviour as for dust in this case too. 
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Next further information on the asymptotics of 9 will be obtained, reach- 
ing the point where the first correction coming from W occurs. Note first 
that (1 + 2r~ 2 W(r cos0)) 1//2 = 1 + 0{r) so that the correction coming from 
W in the second term on the right hand side of the equation for 9 is 0(t~ 2 ). 
The last expression in (fT2l) can be used to replace d/dt{r cos 29) as in the 
case W = and this allows the integral of r sin 29 to be treated. To handle 
the last term in the equation for 9 note that 

r _1 W(r cos 0) cos = -A=* _1 W"(0)-^ (3 cos0 - sin 3 #)+0(r 2 logt) (31) 

6y3 d9 v J 

Furthermore 

4 (3 cos 9- sin 3 9) = 4(3 cos0 - sin 3 9) + Oir 1 ) (32) 
a0 at 

Using Lemma 1 this implies that 9 = —t+0(t~ l logt). The relations f|T5|) and 
ffTB"]) follow. The first place where makes a difference in the asymptotic 
expansions is in the analogue of ffTTl) . It is given by 



9 = -1 



t" 1 sin 2t —t^W'"^) cos 3 t + 0{t' 2 log*) (33) 

2y3 



The results obtained above in the case m — 1 can be generalized to any 
positive m. To get the asymptotic formulae in the general case it suffices 
to replace t by mt everywhere. Theorems are obtained for the potentials 
|m 2 2 + A0 4 and |A(0 2 — 0q) 2 considered in pQ. In the first case the late-time 
asymptotics are proved for all solutions while in the second they are proved 
for all solutions starting sufficiently close to one of the two local minima of 
the potential. 

These results can also be applied to obtain information about the f(R) 
theories of gravity which are equivalent to the Einstein equations coupled 
to a nonlinear scalar field via a conformal transformation. For a discussion 
of this see section 5.6 of [B]. For example consider the theory where the 
Einstein-Hilbert Lagrangian is replaced by R + aR 2 for a negative constant 
a. The potential of the corresponding scalar field is non-negative and has 
a unique minimum at zero. Its second derivative at that point is — 87r/3a. 
Hence the theory just developed can be applied to this case. The interval 
[—01, 0i] and a bound for the value of H at the initial time must be chosen 
appropriately. Similar arguments work for more general choices of f{R). 
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Assume that /(O) = 0, J|(0) = 1 and |^(0) < 0. Then the corresponding 
scalar field has the properties which allow the results of this section to be 
applied to it. Since — > as t — > oo the conformal factor tends to one as 
t — > oo and the physical metric has the same leading order asymptotics as 
the conformally rescaled one. 

It may be remarked in passing that the results of [9] also have applica- 
tions to determining the late-time behaviour of f(R) models with continuing 
acceleration. If there is a positive constant V such that /(— 2V ) = —V 
and /'(— 2Vo) = 1 then the corresponding potential has a minimum at zero 
where it takes the value V . Provided f"(—2V ) lies in the interval (— ^,0) 
then the second derivative of the potential at zero is positive and the results 
of [9J can be applied. Note that the case f(R) = —R 2 discussed in [16J 
(without symmetry assumptions on the solutions) is a borderline one since 
/"(— 2Vo) = — 2^r. In that case the potential is constant. 

Another type of generalization is to fc-essence where the Lagrangian X—V 
of an ordinary nonlinear scalar field, with X = — |V a 0V a 0, is replaced 
by a more general function L(4>,X). The evolution equation for H with 
H 2 = ^-(2XdL/dX - L) becomes H = SnXOL/dX. Note that in a 
spatially homogeneous solution X is non-negative. This Lagrangian will now 
be specialized to the case L(<f>,X) = X — |0 2 — W(4>,X) where W vanishes 
at the origin up to a remainder of third order. Now 

2XdL/dX - L = X + ^<p 2 -2XdW/dX + W. (34) 

It follows that there is a constant C > such that 

C- 1 (<j) 2 + X)<2XdL/dX -L<C(<f) 2 + X) (35) 

for and X small. This means in particular that when and X are small 
H can only be zero when both and X are zero. As a consequence the 
initial data for vanishes so that vanishes everywhere. This case will be 
excluded from consideration and so it can be assumed as in the cases studied 
previously that H > 0. Under the given assumptions on W it follows that 
H is non-increasing. Thus if it starts small it remains small and tends to a 
limit H as t —>■ oo. The equation of motion of is 
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In order to repeat the compactness argument which has been used in other 
cases all that needs to be ensured is that the coefficient of in this equation 
remains bounded away from zero. If the initial value of H is small enough 
this is the case. Passing to the limit in the evolution equation for H it follows 
that 0oo = 0. The equation of motion for 0^ implies that |^(0oo,O) = 0. 
For sufficiently small this means that 0oo = 0. Hence H = 0. 

Passing to polar coordinates in the (0, 0) plane the equation of motion 
becomes: 



sin 9 



\_d_w_ 2X ^w 



dX 



-r cos 9 sin 9 



dX 2 



( dW\ . d 2 W 
- I 1 — — I (3#0) + 



dX 




OW/dX + 2Xd 2 W/dX 2 
1 - dW/dX - 2Xd 2 W/dX 2 

dW/dX + 2Xd 2 W/dX 2 
- dW/dX - 2Xd 2 W/dX 2 



(3Hsin9) + 



d 2 W 
d(j)dX 



d<pdX 



dW 



r sin 9 — r 



+ cos 9 
_ l dW 



dW nv d 2 W 
1 — — 2X 



dX 



OX 2 



(37) 
i 

(38) 



The expression dW / dX+2Xd 2 W / dX 2 , when evaluated at (r cos 9, r 2 sin 2 9/2), 
is 0(r 2 ). From this it easily follows that 9+ 1 is 0{r) and so 9 grows linearly 
with t. The last term in the equation for f is 0(r 3 ), as is the expression 
<]) 2 d 2 W/d<pdX. Thus 



d _! 



\/T27r sin 2 9 + sin 0r 2 ^ + 0(r 

0(j) 



(39) 



The second term on the right hand side of this equation can be written in 
the form |^p-(0, 0) cos 2 9 sin# + 0(r _1 ). It can be concluded as in the case 
of the ordinary nonlinear scalar field that r = 0(t~ r ) and that H has the 
same leading order behaviour as in the case of dust. 



4 Inclusion of matter 

The previous sections were concerned with the Einstein equations coupled to 
a scalar field (possibly nonlinear) without including other fields describing 
ordinary matter. In this section matter satisfying the dominant and strong 
energy conditions will be added. The energy-momentum tensor is the sum 
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of that of a scalar field and that of the other matter. There is no direct 
coupling between the scalar field and the ordinary matter - they interact 
only indirectly via their coupling to the gravitational field. The assumption 
of a spatially flat Friedmann model is maintained. The equation of motion of 
the scalar field, before substituting for H using the Hamiltonian constraint, 
remains unchanged. The Hamiltonian constraint reads 

4-7T . 

# 2 = y(0 2 + 2 + 2U + 2p M ) (40) 

where where pu is the energy density of the ordinary matter. If H is zero 
at any time t\ the initial data for the scalar field vanishes as a consequence 
of the weak energy condition. Then the scalar field vanishes at all times. 
Since the subject of interest here is the effect of a scalar field that case will 
be excluded. It is assumed from now on that H is positive. The evolution 
equation for H is 

1 



H = -An 



4> 2 + Pu + ^trSj, 



(41) 



where trS^ is the trace of the spatial projection of the energy-momentum 
tensor of the ordinary matter. The dominant energy condition implies that 
the right hand side of (|41|) is non-positive so that H(t) is non-increasing. 

Assuming that the solution of the coupled Einstein-scalar-matter equa- 
tions exists globally in the future it can be concluded that H(t) tends to a 
limit H Q as t — > oo. The quantities 0, <fi an d H are all bounded. Hence <fi is 
bounded. The energy-momentum tensor of the ordinary matter is divergence- 
free and it follows that 

dp M /dt + 3H (p M + ^trSw) = (42) 

Thus, by the dominant energy condition, pu is non-increasing and converges 
to some limit p^ > 0. Also tr Sm < 3p M and it follows that H is bounded. 
Differentiating the equation of motion for the scalar field shows that the third 
time derivative of (f) is bounded. Using these facts a compactness argument 
can be carried out as in the case without ordinary matter. Under the same 
conditions it can be concluded that <fi and <f> tend to zero as t — > oo. Moreover 
it follows from (jUJ) that p M + ^ItSm — » 0. Hence trSV/ — * — 3poo as t — >• oo. 
But then pm+^Sm — 2/w If Poo were non-zero then this would contradict 
the strong energy condition. Hence in fact p^ = 0, H = and puif) 
as t — » oo. 
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Consider for simplicity the linear case W = 0. When the equation of 
motion for is written in polar coordinates the resulting equations are 

f = - Vl2^{r 2 + 2p M ) l/2 r sin 2 9 (43) 
9 = -1- v / l2^(r 2 + 2p M ) 1/2 sin#cos# (44) 

Since H tends to zero as t — > oo the argument that 9 grows linearly with 
t also works in the presence of ordinary matter. It is also true that 9 + 1 
remains bounded away from zero. Analysis of the evolution equation for r _1 
shows that r _1 must grow at least linearly and that r = 0(t~ l ). Due to 
its sign the contribution of the energy density of ordinary matter can only 
increase this rate of decay. Unfortunately this does not immediately give a 
decay rate for H. It would seem that to go further it would be necessary 
to know which of the terms r 2 or pu dominates at late times. In the case 
of a perfect fluid with linear equation of state pm = (7 — 1)pm heuristic 
considerations indicate that if 7 > 1 it is self-consistent to require that the 
fluid has a negligible effect at late times while for 2/3 < 7 < 1 (a model which 
satisfies the strong and dominant energy conditions) this is not consistent. 
Rigorous results on the asymptotics in these cases are not available. 



5 Degenerate minima 

Suppose now that the linear scalar field is replaced by a nonlinear one with 
V{4>) = <f) 2n /2n and no additional matter is included. Then the equation for 
H is unchanged and the equation of motion is 

4> + Vl2n~(<j) 2 + ( f ) 2n /n) 1/2 <j) + <j) 2n - 1 = (45) 

As in the case n — 1 the function H is non- increasing and tends to a limit Hq. 
The argument that Hq = goes through without significant change so that 
(j) and tend to zero as t — > 00. For simplicity only the special case n = 2 
will be considered in the rest of this paragraph. Passing to polar coordinates 
leads to the system 

r = r sin 9 cos 9(1 - r 2 cos 2 9) - 3Hr sin 2 9 (46) 
9 = -sin 2 9 -r 2 cos 4 9 -3Hsin9cos9 (47) 

where 

H = Jyr (sin 2 # + ir 2 cos 4 ^ (48) 
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For the solutions being considered here r never vanishes at any time. It 
follows that whenever 9 = kir or (k + l/2)7r for k an integer 9 < 0. Hence 
there are two mutually exclusive possibilities. Either 9{t) —>■ — oo as £ — > oo 
or 0(£) is eventually trapped between 6*i and 6*i + tc/2 where 29i/n is an 
integer. It will now be shown that the second of these cannot occur. To do 
this it is helpful to distinguish between the case where 9i/n is an integer and 
that where #i/7r+ 1/2 is an integer. In the first case 9 is monotone decreasing 
at late times and 9 — > 9\ as £ — > oo. For r small enough 

r>rsm9(^cos6-3Hsm9\ (49) 

This implies that r is eventually positive, contradicting the fact that H — > 
as £ — > oo. Thus this case is ruled out and it can be assumed that 9i/n + l/2 
is an integer. In that case r is eventually monotone decreasing and it can 
be concluded that r = o(l). For any interval of the form [^1,^2] with 9 2 < 
9\ + 7r/2 it is true that at sufficiently late times any solution for which 9 
lies in this interval satisfies a uniform negative upper bound on 9. Thus 
9 — ► 9\ + 7r/2 as £ — » 00 and = o(r). Thus for any e > there is a time £1 
such that for £ > £1 

— sin 2 9 — r 2 cos 4 # — 3iJ sin # cos 9 

< — sin 2 9 — r 2 cos 4 + er | sin 1 1 cos 9 \ 

< - sin 2 9(l- -\ -r 2 cos 2 9 ( cos 2 - - 



2, 

< (50) 

Using the evolution equation for 9 this gives a contradiction. Thus it can 
be concluded that in fact 9(t) — > — 00 as £ — ► 00. It follows that there are 
infinitely many oscillations of the scalar field in this case too. 

There is a remark on p. 242 of [6] that in the case of a nonlinear scalar field 
with potential V(4>) = <j) 2n /2n the solution should be approximated in some 
averaged sense by a perfect fluid with linear equation of state. A statement 
will now be proved which is a concrete realization of this idea. Multiplying 
the equation of motion for a nonlinear scalar field and rearranging gives the 
identity 

= 2 - 6vr0 2 2 -j t {M>+ ^0 2 ) (51) 

In deriving this the equation for H has been used. The information already 
available concerning the solution justifies integrating this relation from some 
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time t to infinity to get 



<pV\(p)ds = J 2 - 6vr0 2 2 ds - <f>(t)<f>(t) (52) 

In the special case of the power-law potential (f>V'(4>) = 2nV((f>). Integrating 
the evolution equation for H and using the fact that H(t) tends to zero as 
t — > oo shows that 

/■oo . 1 

It ^ ds = -^ H ^) ( 53 ) 

Using this in ( 1521) shows that 

roo / pod . \ 

J (j)V'{(j))ds= [J <p 2 ds\ (l + o(l)), t^oo (54) 

It follows that 

J t °°p(s)ds fn- 



{rMds wi; ;i+0(1)) < 55) 

This corresponds to dust in the case n = 1 and to radiation in the case n = 2. 



6 Conclusions 

In this paper rigorous asymptotic expansions for the late-time asymptotics of 
spatially flat homogeneous and isotropic solutions of the Einstein equations 
coupled to a linear massive scalar field are proved. It is shown what similar- 
ities and differences there are when the quadratic potential of this model is 
modified by higher order corrections. The latter results apply directly to give 
information about the dynamics in certain f(R) theories. Basic asymptotics 
are also obtained for fc-essence models where the leading order terms in the 
Lagrangian near the origin agree with those of a linear massive scalar field. 

When ordinary matter such as a perfect fluid is added to the model some 
statements about the asymptotic behaviour are proved but at a certain stage 
in the expansion a competition arises between the scalar field and the other 
matter and the outcome of this is not settled rigorously here. In the case of a 
potential such as V{<p) = 4 /4 with a degenerate minimum it is shown that 
the scalar field is oscillatory in the sense that it has infinitely many zeroes. 
A relation to the radiation fluid is derived. 

As has already been indicated, the asymptotics derived in section 2 can 
be destroyed by the introduction of spatial curvature. When homogeneous 
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solutions with non-vanishing spatial curvature are considered the equation 
of motion for remains the same but the expression for H given by the 
Hamiltonian constraint and the evolution equation for H pick up extra con- 
tributions involving shear and spatial curvature. For Bianchi types I- VIII 
where R < the extra contributions to H and H are positive and negative 
respectively. As in the isotropic and spatially flat case it can be assumed 
without loss of generality that H > while H is non-increasing and tends to 
a limit Ho > 0. A compactness argument as in previous sections shows that 
0, 0, R and the square of the shear a a b(J ab tend to zero as t — > oo. Beyond 
this point it must be expected that the asymptotics of solutions of different 
symmetry classes diverge, as they do in the case without a scalar field. 

To the author's knowledge the only classes of homogeneous solutions 
of the Einstein equations coupled to a massive scalar field whose late-time 
asymptotics have been analysed further than in this paper are the Bianchi 
type I solutions and the isotropic solutions with non-vanishing spatial cur- 
vature. There is a heuristic discussion of the dynamics in these cases in 
PP and [2]. One intuitive consideration for spacetimes containing only a 
massive scalar field similar to that concerning fluids at the end of Sec- 
tion H] follows from the equation dR/dt = —2HR which holds in a ho- 
mogeneous and isotropic model with non-zero curvature and the equation 
d/dt(a a b<J ab ) = —6H(cr ab a ab ) which holds in a Bianchi I spacetime. It sug- 
gests that in the latter case the influence of the shear at late times should 
be negligible while in the former case the curvature should dominate. Indeed 
the claim in p] and [2] is that the late-time behaviour for negative curvature 
resembles that of the Milne model. Of course in the case of positive curvature 
the solution might recollapse. 

In general it appears that it remains to obtain rigorous results for even 
some of the simplest models related to the massive scalar field. It would be 
interesting to see how the oscillations of the scalar field interact with those 
due to the gravitational field which arise in more complicated Bianchi types. 
For these Bianchi types averaging techniques are necessary to analyse the 
case without scalar field. Cf. [TJ, [12], PZTJ- Ideally a unified approach to all 
these problems should be developed. 
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